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The small-scale CMB temperature we observe on the sky is modulated by perturbations that were 
super-horizon at recombination, giving differential focussing and lensing that generate a non-zero 
bispectrum even for single-field inflation where local physics is identical. Understanding this signal 
is important for primordial non-Gaussianity studies and also parameter constraints from the CMB 
lensing bispectrum signal. Because of cancellations individual effects can appear larger or smaller 
than they are in total, so a full analysis may be required to avoid biases. I relate angular scales 
on the sky to physical scales at recombination using the optical equations, and give full-sky results 
for the large-scale adiabatic temperature bispectrum from Ricci focussing (expansion of the ray 
bundle), Weyl lensing (convergence and shear), and temperature redshift modulations of small-scale 
power. The 8N expansion of the beam is described by the constant temperature 3-curvature, and 
gives a nearly-observable version of the consistency relation prediction from single-field inflation. 
I give approximate arguments to quantify the likely importance of dynamical effects, and argue 
that they can be neglected for modulation scales I < 100, which is sufficient for lensing studies and 
also allows robust tests of local primordial non-Gaussianity using only the large-scale modulation 
modes. For accurate numerical results early and late-time ISW effects must be accounted for, though 
I confirm that the late-time non-linear Rees-Sciama contribution is negligible compared to other more 
important complications. The total corresponds to /nl ~ 7 for Planck-like temperature constraints 
and /nl ~ 11 for cosmic- variance limited data to i max = 2000. Temperature lensing bispectrum 
estimates are affected at the 0.2cr level by Ricci focussing, and up to 0.5u with polarization. 



I. INTRODUCTION 



Primordial local non-Gaussianity can provide strong 
constraints on early-universe physics, and the lensing bis- 
pectrum can be used to constrain dark energy models. 
A non-zero CMB bispectrum should be detected by the 
Planck satellite, and it is important to understand the 
physical origins of this signal. 

In principle a non-linear perturbation-theory numeri- 
cal calculation should be able to calculate the expected 
signal (e.g. as first attempted by (lj), but such a cal- 
culation is both conceptually and numerically complex, 
and as yet has not been done convincingly. Fortunately 
however, both lensing and local non-Gaussianity signals 
are dominated by squeezed bispectrum shapes. These 
correspond to large-scale modulations of the small-scale 
power, where the modulating field is correlated to the 
large-scale CMB temperature anisotropics. In the limit 
in which modulating adiabatic curvature perturbations 
are super-horizon at recombination, the local physics in 
each Hubble patch should be identical, and in this case 
analytic arguments can be used to estimate the bispec- 
trum ^ ||. In this paper I re-derive previous results by 
directly relating observed angles on the sky to physical 
scales at recombination. This gives results that are valid 
in a more weakly squeezed approximation, include linear 
transfer and sky curvature effects, and hence are valid 
over a wider range of angular scales so that they can be 
directly applied to data. Including all the relevant effects 
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may be important because of partial cancellations, espe- 
cially for understanding the lensing signal. I also discuss 
a detailed numerical calculation of the dominating lens- 
ing part of the signal, and estimate the scale at which 
complications due to sub-horizon dynamics are likely to 
become important. 

Before getting into a discussion of the bispectrum, 
it is worth briefly remembering the various effects 
that give rise to the observed linear CMB temperature 
anisotropics. Recombination is an approximately equal- 
temperature surface, so observers close to the surface 
would see the same thing everywhere if there were no 
sub-horizon perturbations. However, we see tempera- 
ture anisotropics on all scales because photons from dif- 
ferent parts of the surface get redshifted by differing 
amounts. On scales that are super-horizon at recombina- 
tion the photon rcdshifting effects arc dominated by the 
potential-well Sachs- Wolfe effect, and the late-time Inte- 
grated Sachs- Wolfe Effect (ISW) from passing through 
evolving potential wells along the line of sight. On sub- 
horizon scales overdensities recombine later and are red- 
shifted less and hence appear hot, and at around I ~ 60 
this contribution cancels with the Sachs- Wolfe, so that 
the anisotropics are actually dominated by Doppler sig- 
nals from velocities at last-scattering. As shown in Fig. [y 
there is no scale on which the Sachs- Wolfe limit is accu- 
rate, and only at I <C 60 arc Doppler effects negligible. 
In the region 10 < I < 100 the signal has contributions 
from Doppler, Sachs Wolfe and density perturbations of 
comparable magnitudes, as well as a significant early 
Sachs- Wolfe contribution. For a squeezed-limit result to 
be trustworthy over a useful range of scales it should ac- 
count for these various effects and be not be too sensitive 
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FIG. 1: Power spectrum of the contributions to the total 
CMB temperature anisotropy Ci . The red A 7 contribution is 
the main small-scale contribution from the comoving temper- 
ature perturbations at last-scattering, with $/3 being the net 
large-scale contribution from photons climbing out of poten- 
tial wells (assuming matter domination). Note that the <E>/3 
and A 7 source terms have opposite sign, so their total contri- 
bution to the power spectrum is nearly zero at I ~ 60, where 
the total is then dominated by the Doppler term (green). The 
magenta ISW contributions come from the late-time change 
in the potentials when dark energy becomes important at low 
redshift (solid), and the early contribution (dash-dotted) from 
time-varying potentials soon after recombination as the uni- 
verse became fully matter rather than radiation dominated. 



to the presence of the first corrections from sub-horizon 
dynamical evolution. 

The large-scale temperature anisotropies are gener- 
ated by large-scale perturbations modulating the amount 
of redshifting of a uniform-temperature last-scattering 
surface. The bispectrum arises from the small-scale 
anisotropies on this surface being modulated by the 
large-scale perturbations. If there arc only super-horizon 
modes from single-field inflation present, the physics 
at recombination should be identical in different Hub- 
ble patches, and the small-scale anisotropies would lo- 
cally look statistically equivalent. However when we ob- 
serve the surface these small-scale anisotropies can be 
modulated by large-scale perturbations that change the 
amount of Ricci and Wcyl focussing, Weyl shearing, and 
redshifting that we see in different parts of the sky. Since 
these modulating perturbations also generate the large- 
scale CMB temperature anisotropies, the modulation is 
correlated to the temperature and there is a bispectrum. 

I calculate this signal by integrating along a ray bun- 
dle to relate observed angles to physical scales on the 



last-scattering surface; a concrete calculation of the 'pro- 
jection effects' discussed by Ref. ||. This gives a sim- 
ple derivation of results that are valid on the full sky 
(no flat-sky approximation, which is important for large- 
scale modulations), and results can easily incorporate the 
linear-theory transfer functions to include all the relevant 
terms in the modulation correlation function. Previous 
work has cither focussed on the CMB lensing bispectrum 
alone, or has done a full approximate flat-sky calculation 
neglecting the part of the lensing bispectrum that is not 
due to super-horizon modes. However this separation is 
not quantitatively justified, since the lensing bispectrum 
has significant contributions from early and late times 
(from early and late ISW effects), as well as large scale 
modes and other sources. The sources do not separate 
cleanly in either redshift or wavenumber, so it is best to 
do a full self-consistent analysis. 

Previous work has emphasised that there are partial 
cancellations between Ricci focussing and the lensing 
(non-ISW) Weyl focussing terms [^|, [|, 1| , and a full anal- 
ysis is therefore potentially important for lensing bispec- 
trum studies (as well for a consistent bispectrum calcu- 
lation for subtraction when looking for primordial non- 
Gaussianity) . This cancellation is very similar to that 
in the linear CMB anisotropics, where the total large- 
scale anisotropy is smaller than individual terms (see Ap- 
pendix |a|). 

The Ricci focussing modulation accounts for the well- 
known 'consistency relation' result for the bispectrum in 
single- field inflation || calculated here in the context 
of the directly-observable CMB, which has acoustic oscil- 
lations, the signal is rather larger than it appears in the 
primordial curvature perturbations (which have a very 
smooth spectrum so the signal is completely negligible). 

In practice the (Weyl) lensing signal is often most 
easily calculated by reconstructing the lensing potential 
using quadratic maximum likelihood estimators [g| (see 
Ref. [§ for a review), and then correlating it with the 
temperature to get the bispectrum: in the squeezed limit 
the temperaturc-lcnsing reconstruction power spectrum 
contains the same information as the bispectrum. By 
using a lensing potential, shear and convergence effects 
are treated self-consistently in one go since they are re- 
lated by gradients at linear order for scalar perturbations. 
For this reason I will group the lensing effects separately 
from the others, and quantify the importance of the par- 
tial cancellations for the full result. This allows biases 
in a pure-lensing analysis to be assessed and corrected, 
and to quantify whether the other terms are significant. 
My classification is therefore closer to Ref. (J] — who 
exclude Weyl lensing altogether — than Ref. [9, who 
include early Weyl lensing in their combined result but 
don't include the larger late-time signal. The full result 
can easily be obtained by adding the terms together. I 
consider the full result including polarization in the final 
section. 

The approach followed here is similar to the analysis 
of redshift-space densities in Refs. [10 . Ref. |It| also 
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recently developed general results for distortions to stan- 
dard rulers in redshift space. The main difference for the 
CMB is that we are now concerned with a single constant 
temperature recombination surface (rather than constant 
redshift surfaces), and that the CMB is at much higher 
redshift so the line-of-sight impact of super-horizon per- 
turbations is relatively more important. 



A. Jacobi map between physical scales and 
observed angles 

For single-field inflationary fluctuations, there is only 
one clock (e.g. set by the temperature), so locally physics 
is the same at recombination in all patches that differ 
only by super-horizon modes (up to gradients that we 
consider separately later) j^, |^, [l2|. However, the small- 
scale physical scales at recombination arc observed at 
different angular sizes and temperatures because of the 
large-scale perturbations. Differences in redshifting from 
the uniform temperature surface in different directions 
gives rise to the linear CMB anisotropics (see Appendix |A| 
for a brief review) . Similarly the same physical sizes ap- 
pear at different angular sizes in different directions, due 
to bending of the photon path and differences in expan- 
sion along the line of sight. The path deflections also give 
rise to shear distortions. 

I calculate the anisotropy in observed angular sizes and 
shear by using the optical equations for the propagation 
of a ray bundle through the perturbed universe. For a 
pair of rays separated by an infinitesimal angular coor- 
dinate 89 j about an observation direction ii at 4, the 
physical infinitesimal separation vector orthogonal to the 
direction of propagation at affinc parameter A along 
the ray is given by 

^{X)=Vu{X)56j, (1) 

where Vjj is called the Jacobi map (see Ref. jl3| for de- 
tailed review in similar notation). In general T>jj can 
be found by integrating the optical equations along the 
line of sight |Q. The trace of Vjj determines the (lin- 
earized) magnification and the symmetric trace-free part 
determines the shear. Scalar perturbations generate no 
rotation at first order so the matrix is symmetric. 

For simplicity I restrict to a flat, almost-Friedmann- 
Robertson- Walker (FRW) model with scalar perturba- 
tions, working to first order in the conformal-Newtonian 
gauge (CNG) with metric 

ds 2 = a 2 (ry)[(l + 2^)dr/ 2 - (1 - 2$)*yda:W]. (2) 



perturbed universe 1 (see e.g. Rcfs. [^0|, |l~5| |) 
V(h, rj)/2 = x(n, r])a(v) [1 + $a - $ - k + n • v A ] ■ (3) 

Here the lensing convergence k is determined solely by 
the Weyl focussing along the line of sight, va is the 
Newtonian-gauge (peculiar) velocity of the observer at 
A, and the local potential &a is included so that the ob- 
server's scale factor is defined to be unity, a(0)(l — &a) = 
1. The potential at the source $ looks like a local term, 
but is just the result of integrating a total derivative part 
of the Ricci tensor along the ray: it accounts for the line- 
of-sight Ricci focussing that is not accounted for in the 
convergence term n generated by Weyl focussing. 

The Riemann tensor only depends on second deriva- 
tives of 4> and so only these gradients are observable: 
a constant <J> cancels in Eq. (||) so only differences in po- 
tentials are observable. A constant V 2 $ gives a physical 
monopolc contribution to T> but is equivalent to the pcr- 
turbative effect of a non-flat FRW geometry p5[ , so I 
shall not consider the monopolc further and drop $a- 
As is usual in CMB lensing studies we can then define a 
lensing potential for a source at radial distance x* 

i, = -2 dx^^Vw(xn,VA-x) (4) 
Jo XX* 

where the Weyl potential = (<E> + \&) /2 determines 
the scalar Weyl (conformally invariant) part of the Rie- 
mann tensor, so the deflection angle is V itp and the lens- 
ing convergence and shear are given by k = — Vi\7 I ip/2, 
7jj = V/jVj)^. I choose to define 'lensing' to be the 
Weyl effects, not including the Ricci focussing term that 
only depends on differences in potentials. 

For modulation scales that are large compared to the 
thickness of the last-scattering surface it is a good ap- 
proximation to treat the recombination visibility as a 
delta function at time 77, , a distance x* away in an FRW 
universe. Recombination happens when the universe 
cools enough so that hydrogen can recombine, so in this 
approximation last-scattering is an equal-temperature 
surface. We can then consider how small-scale sub- 
horizon perturbations on this approximately equal tem- 
perature surface are observed, allowing for the fact that 
there are large-scale perturbations that do not change 
the local physics at the source — see Fig. |[ If 77* is the 
background (conformal) time of recombination, in the 
perturbed universe recombination will happen at a per- 
turbed conformal time 77 = 77* + Srj (on sub- horizon scales 
this is the physical effect that overdensities need longer to 
expand to the recombination temperature): for an equal 
temperature surface with fixed photon energy density we 
need pL8r} = — 8p 1 and hence 5a/ a = T-LS-q = where 
A 7 is the Newtonian-gauge fractional photon density per- 
turbation. In terms of the background scale factor a* and 



The trace of the Jacobi map T> depends on the Ricci 
tensor and determines the observed angular size; it gives 
a generalization of the angular diameter distance to a 



1 This assumes the source plane is ray-orthogonal, however the 
correction is next order so we can neglect it. 
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FIG. 2: Schematic equivalent descriptions of an observation of fixed physical-size small-scale perturbations on the equal- 
temperature recombination surface, in the presence of large perturbations that were super-horizon at recombination. The 
perturbation on the right side of each figure appears to have slightly smaller size on the sky. Ricci (de)focussing directly 
changes the cross-section of the ray bundle due to the expansion of spacetime along the beam. Weyl lensing bends the path of 
the beam, causing Weyl focussing (and shear) by differential deflection. Angular anisotropy in the Ricci focussing is described 
gauge-invariantly by the constant temperature 3-curvature. In the left figure, a non-comoving gauge with constant temperature 
slices, Ricci (de)focussing is what changes the observed angular size of the perturbations. In the comoving Newtonian gauge on 
the right figure, the effect can be described as the overdensity having recombination happening later, so fixed physical scales 
have smaller comoving size and appear smaller in the background, combined with a perturbation to the scale factor at the 
unperturbed background comoving distance described by the potential. The Ricci focussing effect that makes the patch look 
smaller is then partly compensated by Weyl focussing due to differential deflection of the beams in the potential gradients along 
the line of sight; this cancellation is not shown to more clearly illustrate the Ricci effect. The right patch also looks colder 
because there is more redshifting from the photons climbing out of the deeper potential well than the partly-compensating 
blue-shifting from the beam contraction (see Appendix Q note sub- horizon overdensities would look hotter). 



radial distance x*i the physical scales on an equal tem- 
perature surface are then related to angles to the sky by 



T>/2 = x* a * 



X* 4 



n - v A 



(5) 



Recall that recombination happens at a time 77* ~ 
270Mpc and the distance to last scattering in a stan- 
dard ACDM model is x* ~ 14Gpc, so ?7*/x* ~ 0(0.01). 
The 0(A 7 r]*/x*) term in the radial displacement Sx/x* 
is therefore much smaller than the others and can be 
neglected, and the time delay contribution is also very 
small due to cancellations The effect of lens- 

ing convergence dominates because although the poten- 
tials are roughly constant after recombination the sig- 
nal is enhanced (even incoherently) by the large number 
0(x*/r]*) °f perturbations along the line of sight; the 
convergence signal is also blue because of the two an- 
gular gradients. However, there is only a correlation of 
the convergence with the CMB temperature (and hence 
a bispectrum) from very-large scale lenses, lenses close 
to the last scattering surface, and the integrated Sachs- 
Wolfe effect (see later section for detailed discussion). 
The term — <T> is of order WSri and significantly smaller, 
but highly correlated to the temperature and potentially 
marginally detectable in the bispectrum (and hence also 
potentially leading to biases if neglected) : one of the aims 



of this paper is to quantify more carefully the importance 
of this term. The local Doppler aberration term n • va 
is relatively large 0(1O -3 ), but only contributes a dipole 
modulation and is not usually considered to be part of 
the cosmological signal since it depends on the local pe- 
culiar velocity. 

In the background the angular diameter distance 
T>/2 = x* a * is the Ricci focussing result for a flat FRW 
universe (Weyl is zero in the background). A 7 /4 is the 
perturbation to this background result due to differences 
in time coordinate (hence perturbations in the expan- 
sion). The $ term is the perturbation in the Ricci fo- 
cussing to the background recombination time, and the 
sum A 7 /4 — <E> gives the total perturbation in observed 
size due to Ricci focussing. The decomposition is not 
gauge- invariant, but the sum £ 7 = A 7 /4 — $ is the 
gauge-invariant 3-curvature on constant photon density 
hypcrsurfaces. In a sense the Ricci focussing pertur- 
bation is the ray bundle area analogue of the SN for- 
mula for the comoving curvature perturbation from in- 
flation: differences in expansion between different Hub- 
ble patches in inflation are seen as varying expansions of 
the ray bundle between the observer and different points 
on the last-scattering surface. An overdensity at recom- 
bination has larger local scale factor (described by the 
3-curvature) than the reference point, so the beam cross- 
section shrinks compared to average as the ray bundle 
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FIG. 3: The CMB temperature cross-correlation power spec- 
tra relevant for contributions to the bispectrum (dashed lines 
are negative). The blue curve is the 'standard' result for 
lensing convergence and shear, magenta shows the contribu- 
tion from Ricci focussing determined by the super-horizon 
curvature perturbation £q without evolution (the correlation 
changes sign at I ~ 60 where Sachs- Wolfe and comoving over- 
density terms cancel in the temperature). At / > 100 sub- 
horizon effects may become comparable, with terms estimated 
for large-scale flows (red dashed) and perturbations in sound- 
horizon (cyan). The solid red curve shows an estimate of the 
total effective magnification correlation including flows and 
Ricci focussing. The total effective shear is only modified by 
sub-horizon dynamics, and is estimated by the dash-dot blue 
line. The sub-horizon effects are unlikely to be accurate, but 
give an indication of the scale on which they might become 
relevant. The shaded band shows an indicator of ±1<t sen- 
sitivity to systematic errors rel ative to the cosmic variance 
limit from lensing, given by ± v /C i TT C i KK + (Cf K ) 2 /I. 



leaves the perturbation, leading to fixed physical sizes 
appearing smaller (see Fig. |J). 

The result for V/ 2 is valid on all scales, but the approx- 
imation of constant physical scales on the source plane 
may be modified by dynamics. For the limiting case of 
super-horizon adiabatic modes where the physics is iden- 



tical in single-field inflation, 



$ = £, where C is the 



conserved 3-curvature perturbation on equal (total) den- 
sity hypersurfaces (in matter domination £ w — 53>/3). 
Note that overdensitics have positive k, positive A 7 and 
negative <f>, so lensing convergence tends to cancel with 
the other terms: an overdensity has Ricci defocussing 
compared to average, but also acts as a magnifying lens, 
so the effects partly cancel. For a scale-invariant spec- 
trum in matter domination and no ISW, in the flat- 
sky approximation 5/6 of the effective magnification- 



temperature correlation is cancelled [|2|, |5J , so it is impor- 
tant to consistently model both effects unless each one is 
separately negligible compared to the ISW signal. Note 
that a constant super-horizon A 7 perturbation appears 
to change T>, however this is correct because it would be 
equivalent to observing a larger CMB background tem- 
perature, which for the same physics requires less time 
since recombination and hence corresponds to a smaller 
comoving radial distance to recombination. 

The shear is significantly simpler since it vanishes in 
the background, and is given by the symmetric-trace free 
part of the Jacobi map T^uj) = cl*X*1Ui where yu is 
given in terms of the lensing potential as described above 
(see e.g. |jl3|). The shear is given solely by Weyl effects 
(lensing) , so there are no complications with partial can- 
cellations. 

The shear and convergence are angular derivatives on 
the sky, and the lensing can be thought of as re-mapping 
points by a deflection angle (which is a conveniently sim- 
ple effect to simulate). However, the Ricci focussing is 
not an angular derivative, and cannot be modelled using 
deflection angles: for non-trivial curvature of the recom- 
bination surface there is no local-isotropy preserving map 
onto the observed sphere. The spherical projection is in- 
trinsically 3-dimensional, with the Ricci term being more 
like a radial deflection than an angular one. 



B. Sub- horizon dynamics 

Leading effects of dynamics may be important, because 
velocities come to dominate the CMB power spectrum 
at I ~ 60 (see Fig. [I]) and so only the very lowest multi- 
poles can safely be modelled as due to having only static 
super-horizon sources at recombination. Dynamics will 
modify the small physical scales observed, and hence lead 
to additional anisotropy in the distribution of the small- 
scale perturbations on the last-scattering surface. Here 
I give some simple arguments to quantify these effects, 
and hence determine a range of scales over which a simple 
squeezed analysis is likely to be robust. 

Note that from the photon density perturbation evo- 
lution equation 



3> 



-J V.v 7 d,/ 



Co - r / V • v 7 d??' 



(0) 



for adiabatic perturbations. We can therefore isolate 
the dynamical velocity compression at recombination by 
defining 



SN V = HSr] v 



(7) 



which gives the fractional change in expansion due to 
the large-scale flows in the fluid compressing the small- 
scale perturbations, with 8rj v being the corresponding 
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change in conformal time until recombination. Note 
that it is small on super-horizon scales since 5N V ~ 
0((k/H) 2 Co). For scalar perturbations the flow is po- 
tential v = Vi> 7 , so there is a deflection angle model 
a v = — J dr]v±/x* = — Vn / drjv^/xt for the dynamic 
effect of the bulk velocity-induced motion of small-scale 
perturbations in the plane of the sky, assuming the 
small perturbations co-move with the large-scale radia- 
tion flow. This gives rise to an effective shear and conver- 
gence compared to purely static evolution. However, the 
3-divergence in flow increases the local temperature; this 
is removed by additional expansion, so the total effect on 
the physical scale at the equal-temperature recombina- 
tion surface is only the net effect. Defining the angular 
velocity distortion potential 



'* drjv^ 



V~ 2 SN V 



X.* 



X.* 



(8) 



we can then calculate angular convergence k„ = 
—ViV I ip v /2 and shear rfj = V(jVj}^ with respect 
to primordial isotropic initial conditions, in exactly the 
same way as we do for lensing using the lensing po- 
tential. The total effective fractional magnification and 
shear with respect to isotropy, including Ricci, Weyl and 
angular velocity dynamics, is then roughly 



-eff 



Co 



(9) 



In addition to the bulk movement of the small scale- 
perturbations, there will also be a change in their corre- 
lation length, determined roughly by the sound horizon. 
The fractional change in sound horizon due to recombi- 
nation being delayed in a hotspot gives a scale magnifi- 
cation ~ Srjv/r)* which reduces the effect of the velocity 
convergence but does not alter the shear. This is only 
a rough order of magnitude as there are of course many 
other effects, for example due to ionization fraction per- 
turbations on the constant temperature slices (c.f. [p7|). 
There will be additional velocity effects from the radial 
components of the velocity, however these are more like 
power modulations similar to /nl, and are therefore not 
potentially enhanced by d In Ci/d In I due to the acoustic 
structure, and should be safely negligible at least as a 
contaminant of the large-scale lensing cross-correlation. 
There will also be an /NL-like signal from modulation in 
power due to changes in the Silk scale and recombination 
thickness. 

Figure |3| shows a numerical calculation of angular cor- 
relation power spectrum of the various contributions with 
the CMB temperature using CAMB |^8| . The effective 
magnification correlation is enhanced on large scales be- 
cause there is a near-cancellation between convergence 
from large-scale lenses and the Ricci focussing (o term; 
the large-scale modes contribute to the correlation with 
opposite sign to the ISW signal, so the signal is enhanced 
when most of the large-scale modes are cancelled by the 
Co modulation. 

Note that temperature anisotropy directly probes the 
velocity via the Doppler effect which is O(krj^), but flow- 



generated convergence and shear depend on the effect 
of gradients of the velocities over time and hence are 
0((/c?/*) 2 ). It is therefore required to model the evolution 
and Doppler contributions to the large-scale temperature 
sources on significantly larger scales than is required for 
dynamical convergence and shear, or any other effect that 
depends on changes in time or temperature due to fluid 
compression. Furthermore the Doppler contribution to 
AT is uncorrelated to scalar densities, so the correlation 
of AT with dynamical compressions is small around I ~ 
60. Only on smaller scales, roughly I > 100 does AT 
start to become dominated by the flow-induced density 
perturbations, and at that point a full analysis would be 
required to get an accurate result for the correlation. 

These arguments and the numerical estimates indicate 
that a squeezed analysis neglecting second-order dynam- 
ics is likely to be accurate at I < 100. The 0((kr]) 2 ) scal- 
ing of sub-horizon effects should ensure that the scaling of 
the full result is similar to the approximations here, and 
hence that this conclusion is fairly robust. Since I < 100 
contains most of the ISW signal this is likely to be a good 
data cut for parameter studies that want to be confident 
of avoiding issues with sub-horizon dynamics. For pri- 
mordial local non-Gaussianity the signal goes to zero at 
a modulation scale I « 60 where G l changes sign, so a 
restriction to I < 60 should ensure dynamical effects are 
small while losing little signal (for Planck ~ 90% of the 
temperature signal is at I < 60, see e.g. Ref. jl9[ Fig. 11 
]). A cut at I < 100 is also sufficient for corrections due 
to finite thickness of the last scattering surface not to be 
large. 



II. THE CMB BISPECTRUM 

If we observe a small patch of the CMB, the statis- 
tics of the temperature fluctuations may vary as a func- 
tion of position on the sky due to large-scale magnifi- 
cation, shear, and temperature modulation. Correlation 
of these modulations with the CMB temperature gives a 
bispectrum. If the (non-linear) modulated temperature 
is T m , for l\ < lii h an d neglecting tiny non- linear effects 
on the large-scale mode (the linear short-leg approxima- 
tion), the bispectrum from a set of Gaussian modulation 
fields Xi is given by 



(T imi T; 2m2 T; 3 „; 3 ) ~ C (i 



TXi 



SX. 



(10) 

See Ref. J20| for a derivation of this non-perturbative 
weakly-squeezed version of a well-known result: the bis- 
pectrum depends on the large-scale correlation of the 
modulation with the temperature, and the response of 
the small-scale power to changes in the large-scale mod- 
ulation modes. Only the lensing effect is significantly 
large, so the non-perturbative average is only required 
over temperature and lensing potential modes. This gives 
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results that (approximately) involve the lensed small- 
scale power spectra rather than the unlensed ones. 



A. Weyl lensing bispectrum 

The largest term is the lensing bispectrum, which is 
well understood Jl9|, ^T|-p4|. The temperature bispec- 
trum signal is limited for detection by cosmic variance 
at about 5a, and should soon be detected by Planck at 
around 4<r. From Ref. Jl9) the reduced bispectrum is 
given to good accuracy on all scales by 

h Ms = \ [(h(h + 1) + h(h + 1) - h(h + 1)] cl 4 'C l2 

+ perms (11) 

where ip is the lensing potential. Note that this is non- 
perturbatively accurate in the small-scale modes, since 
the bispectrum involves the lensed small-scale power 
spectra - an effect that would be missing in any complete 
purely second-order-only analysis. In the flat-sky ultra- 
squeezed limit li <C h,h this reduces up to 0((li/l) 2 ) 
corrections to 



J hhh 



1 djPCt) 
I 2 d\nl 



cos2(f> h i 



dQ 
dhd 



(12) 



where 1 = (I2 — h)/2 and 0;^ is the angle between 
li and 1 (the full sky version of this limit is given in 
Ref. |2f| appendix B]). The first term describes the 
isotropic squeezed bispectrum from magnification, and 
the second term the anisotropic quadrupolar bispectrum 
from shear. In the squeezed limit the two terms are effec- 
tively orthogonal. If neglected the temperature lensing 
bispectrum on its own is known to bias local primordial 
/nl estimators by /nl ~ 6 for Planck (/nl ~ 9 in the 
noise- free limit) |24|, ^7j, but the distinctive shape 
makes lensing easily distinguishable. 

Dynamical potential flow of small-scale perturbations 
gives an effective convergence and shear as described in 
Sec. IB, and the form of the induced bispectrum from 



these effects will be identical to lensing. The local n • va 
Doppler aberration in Eq. (j^) also contributes the same 
way as dipolar convergence (there is no shear dipolc) [|13| , 
and hence gives a bispectrum of the form of Eq. ( |ll] ). 
Since v,i also determines the direction of the CMB tem- 
perature dipolc, in a sense there is a large local bu 2 i 3 bis- 
pectrum that should be detectable by Planck [^8j , though 
this is not usually regarded as cosmological as it depends 
on the local velocity. 



in Eq. (||) is equivalent to the effect of distance per- 
turbations when looking through an unperturbed uni- 
verse to a last-scattering surface at equal time, with 
2?/2a* = x*(l + Co)- In this conformal picture the change 
in angular scale is coming purely from changes to the ef- 
fective distance away 2 , and the reduced bispectrum from 
the Co distance modulation is then easily calculated to be 



°hhh 



1 



2 dlnx* 



Cu + Cu 



(13) 



For a scale-invariant C\ spectrum with sharp recombi- 
nation Ci is independent of \*-, so there is no visible 
scale modulation and hence no bispectrum, as expected; 
however in reality C\ has significant acoustic oscilla- 
tions. For the ultra-squeezed limit with l\ <C 200, taking 
1(1 + 1)C; oc f(l/x*) this reduces approximately to 



J hhh 



-C 



TCo 



1 d 
I 2 dlnZ 



(14) 



which as expected has the same form as the isotropic con- 
vergence term in Eq. (|l2| ) . Note that cf^° goes through 
zero at l\ ~ 60, so it should not be approximated with 
the h < 60 Sachs- Wolfe result. 

Eq. @ is the CMB equivalent of the - V{kx)^^ 
'consistency relation' bispectrum generated by single- 
field inflation @, Note however that it is not correct 
to evaluate an effective /nl = 0(n s — 1) in the primor- 
dial power spectrum and then use that in the result for 
physical small-scale power modulations (Eq. ( |l6| ) below) , 
the reason being that d/dlnfc does not commute with 
the CMB transfer functions. Since the transfer functions 
arc highly oscillatory the signal from Eq. ([l4]) is much 
larger (though still small). If there were no other signals 
it would be (un) detectable 3 at ~ 0.4<r (using l\ < 100), 
and projects onto /nl as a bias of /nl ~ 1. If neglected 
it also contributes to the lensing signal, corresponding 
to a bias of around 5% in lensing bispectrum amplitude, 
corresponding to a bias of about 0.2o\ However it is also 
easily calculable so these small biases can easily be mod- 
elled or subtracted. 

The correction to the squeezed trispectrum is safely 
small compared to lensing, equivalent to only 1 — 2% on 
the k power spectrum on the very lowest multipoles — 
well below cosmic variance. 



C. Redshift and power modulation bispectrum 

In the single-surface approximation, the perturbed re- 
combination surface is equal temperature. However, we 



B. Ricci focussing bispectrum 

The Ricci focussing term gives an effective observed 
magnification of the observed angular scale (but no 
shear). In a conformal picture the effect of £ ~ A 7 /4 — $ 



2 Note that in the equivalent conformal picture overdensities are 
further away, which is opposite to the expanding picture where 
overdensities recombinc later and are actually slightly closer in 
the absence of large time delays. 

3 I calculate biases, correlations and errors using a minor modifi- 
cation of CAMB's public full-sky bispectrum code |l9|. 
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see temperature anisotropies because the photons get dif- 
ferent total amounts of redshifting along different lines 
of sight. When we observe small-scale anisotropies, the 
observed photon temperatures are also modulated by 
the large-scale variations in redshifting between differ- 
ent points on the last-scattering surface. Note that this 
is exactly the same effect that gives rise to the linear 
large-scale anisotropies (Appendix , so the small-scale 
anisotropies are just modulated by the large-scale 1+AT. 
Allowing for the fact that the small-scale power is also 
lensed this gives a reduced bispectrum 4 



Ci 1 ( Ci 2 + Ci, 



(15) 



Local primordial non-Gaussianity (for example from 
multi-field inflation) also generates local physical mod- 
ulations in the small-scale perturbations, conventionally 
parameterized as ( = [1 + (3/NL/5)Cg]Cg- The CMB bis- 
pectrum from this signal is well known (see e.g. [^9|-j3lf ), 
and the full result can be easily calculated using CAMB. 
To compare to the other weakly squeezed signals, we can 
take the limit where l\ corresponds to a scale much larger 
than the thickness of the last-scattering surface (e.g. the 
ultra-squeezed limit). The bispectrum then reduces to 
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D hhh 



(16) 



where Co is the primordial curvature perturbation evalu- 
ated at the position on the last-scattering surface (a nu- 
merical calculation of Cf° is included in Fig. ||). This 
approximation is not accurate at l\ > 60 due to finite 
recombination thickness, but is nonetheless ~ 90% cor- 
related with the full result and is accurate for low l\. 
Note that it also trivially includes the effect of lensing 
on the primordial non-Gaussianity because the modula- 
tion affects the lensed small-scale power spectrum, see 
Ref. [||. 

On very large scales, neglecting ISW, AT ~ — Co/5 so 
in the ultra-squeezed limit bi 1 i 2 i 3 ~ — 12/nlC'; 1 C;; hence 
the always- present signal of Eq. ( |l5|) is comparable to 
/nl ~ — 1 /6 and therefore negligible compared to cosmic 
variance limits of /nl ^2 @ . The shape of Eq. ( |l5|) is 
actually different to local non-Gaussianity because C z 
has opposite sign to Ci x for l\ > 60, but remains even 
less detectable than the Ricci focussing effect because it 
is not enhanced by oscillations due to dlnCz/dlnZ fac- 
tors (which have amplitude ~ 6). Note however that the 
local dipole part of the temperature modulation gives an 
easily detectable dipolar power modulation signal p8|, 



4 Note that Ci l here includes the ISW effect: viewing small-scale 
perturbations through a large evolving potential well will give a 
large-scale modulation of the observed redshift due to all relevant 
large-scale effects including ISW (the authors of Ref. [bl agree 
with this contrary to a statement in the text). 





bias on /nl 


Data used 


°7nl 


Weyl 


Ricci 


Redshift 


Total 


T 


4.3 


9.5 


1.5 


-0.22 


10.7 


Planck T 


5.9 


6.4 


1.0 


-0.22 


7.1 


T (h < 60) 


4.6 


10.6 


1.7 


-0.25 


12.0 


Planck T (h < 60) 


6.2 


7.0 


1.1 


-0.25 


7.9 


T+E 


2.1 


2.6 


1.1 


-0.05 


3.7 


Planck T+E 


5.2 


4.3 


1.0 


-0.15 


5.2 



TABLE I: Individual and total biases on primordial local- 
model non-Gaussianity parameterized by /nl for CMB tem- 
perature and _E-polarization data with Planck-like noise (as- 
suming isotropic coverage over the full sky with sensitivity 
AT = AQ/2 = AL/2 = 50MKarcmin [ A/ ; T = iVf/4 = 
2 x 10~VK 2 ] and a beam FWHM of 7 arcmin) or cosmic- 
variance limited data with l max = 2000. Results are assuming 
that non-/NL contributions are only significant at l\ < 300 
and negligible dynamical effects; the h < 60 results are fil- 
tered to only use large scale modulations and are therefore 
immune to small-scale modulation effects. The bias is the 
systematic error on /nl if the given contribution is neglected, 
which can be compared to cr/ NL which is the Fisher error es- 
timate (including lensing signal variance). 



which must be carefully subtracted to get at a primor- 
dial tnl trispectrum pa ] (along with the dipole aberra- 
tion discussed in the previous subsection). 



D. Full squeezed bispectrum 

Even if there is no primordial non-Gaussianity, there 
is still a detectable CMB squeezed bispectrum due to 
the effect of perturbations on light propagation. This 
signal is all that is expected in single-field inflation mod- 
els, and must be modelled to constrain primordial non- 
Gaussianity. The full signal is given by summing the 
three terms calculated above. The resulting projec- 
tions onto primordial /nl are shown in Table | The 
Ricci focussing term also affects the Weyl lensing bis- 
pectrum, corresponding to a bias of 0.04 for both Planck 
(fi G ns = 0.26) and cosmic variance (<7i ens = 0.19): a bias 
of ~ 0.2o\ 

Including all the terms the projection on /nl increases 
from /nl ^6.4 for lensing only to a combined signal of 
/nl ~ 7.1 for a Planck- like temperature measurement 
with the same assumptions as Ref. |lj|. The increase 
is due to the cancellation between the large-scale Ricci 
focussing and Weyl convergence, which acts to increase 
the total magnification-like signal compared to that cal- 
culated from lensing deflection angles alone. This is be- 
cause the large-scale modes contribute to the tempera- 
ture correlation with opposite sign to ISW. The correc- 
tion is however small compared to the expected error bar 
of cr/ NL ~ 6. The projection onto the lensing bispectrum 
is somewhat larger because Ricci focussing looks like con- 
vergence, and the potential bias is around 0.2rr. 
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FIG. 4: Contributions to the CMB lensing deflection angle- 
temperature cross-correlation power spectrum, defining the 
lensing potential using a recombination visibility weight func- 
tion. Dashed lines are negative and the shaded band gives 
an estimate of cosmic variance sensitivity as in Fig. |5| The 
results are not directly observable on scales I 3> 100 where 
recombination surface effects become important (and hence 
acoustic oscillations), but indicate the order of magnitude of 
the complicating effects. At I < 100 where the signal is un- 
der good control it is dominated by early and late ISW, with 
a negative contribution from very large-scale modes at the 
lowest I and a small Doppler contribution. The non-linear 
Rees-Sciama signal is always small compared to other effects. 
The CAMB Jan 2012 result using a delta-function visibility ap- 
proximation and cutting sources very close to recombination 
agrees well with the black curve for / < 100. Note that this 
plot differs by a + 1) scaling factor compared to other 
plots in this paper. 

This is assuming that there are no dynamical sub- 
horizon effects that project onto /nl; a conservative 
analysis can restrict to l\ < 60 where the squeezed 
results should be robust, and dynamical and finite 
recombination-thickness lensing effects are small. For 
Planck this restriction only increases the error bar by 
~ 5% (conversely any dynamical effects would have to 
be very large to significantly bias a full analysis, since 
only ~ 10% of the temperature signal is from higher li). 



III. TEMPERATURE-LENSING CORRELATION 

In the absence of primordial non-Gaussianity, the dom- 
inant contribution to the CMB bispectrum is CMB lens- 
ing. It is therefore important to model it accurately. 
Eq. ([n]) is an accurate result for the lensing bispectrum, 
but requires a calculation of Cj^' , the correlation between 
the temperature and the lensing potential. Here I inves- 
tigate this signal in more detail to quantify the relevant 
physical contributions, assess where the calculation is re- 



liable, and whether non-linear Rees-Sciama contributions 
are important. 

The following contributions may be relevant for the 
temperature lensing-potential correlation: 

• Large-scale modes that span recombination and 
also act as lenses 

• The late Integrated Sachs Wolfe effect (late ISW) 
at low redshift from decaying potentials 

• The early Integrated Sachs Wolfe effect (early ISW) 
due to the transition from radiation to matter dom- 
ination, and decaying modes 

• Lenses close to last-scattering being correlated to 
density perturbations that have infall giving a 
Doppler signal in the CMB 

• Doppler signal from scattering at reionization 

• Lenses at last-scattering that directly correlate per- 
turbations to lensing at the recombination surface 

• Non-linear Rees-Sciama signal at low redshift from 
non-linear gravitational clustering 

• Non-linear SZ signal from scattering in clusters 

• Correlations due to foreground contaminants 

The signal of most interest for parameter studies is the 
linear late ISW since it can be used as a constraint 
on dark energy models j32| - the CMB lensing po- 
tential happens to correlate at the 90% level with the 
ISW (3j|. However the total temperature correlation 
drops off rapidly with scale as the ISW component of 
the signal tends to cancel between multiple perturba- 
tions along the line of sight. Except for the last three 
non-linear and foreground signals, all of the other ef- 
fects are self-consistently calculated in linear theory by 
using a full standard linc-of-sight Boltzmann code, and 
are included CAMB's standard calculation 5 . The thermal 
SZ signal is frequency dependent, but depending on fre- 
quency potentially important at I ~ 100 |34(| : I shall not 
consider this signal further here since it is zero at the 
SZ null and in principle is isolable by its frequency de- 
pendence, as are other potentially significant foreground 
signals like the correlation with the cosmic infrared back- 
ground (CIB) J35], |36| (the CIB signal is also small at 
frequencies usually used for cosmological analysis). 

The late-time ISW effect can also be calculated accu- 
rately on intermediate and small scales using the Limber 



5 CAMB Jan 2012 uses a single source-plane approximation and 
cuts sources very close to last-scattering; the result agrees well 
with the total result shown here at / < 100 but does not in- 
clude most of the recombination surface effects that cannot be 
modelled consistently anyway. 
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approximation result |32| 



rc, 



X* 



X 9 



drj 



V*{k = l/x,v), (17) 



where is the power spectrum of the potentials that 
can easily be calculated from the standard matter power 
spectrum. The approximation of using <9, ) 'P$ is not 
strictly valid for non-linear evolution, but has been tested 
on simulations and appears to be adequate j3?j , in which 
case the same result can also be used to calculate the 
full late ISW+Rccs-Sciama result simply by using the 
non-linear power spectrum. I use the Halofit [|38| non- 
linear matter power spectrum model to estimate the dif- 
ference in the total matter power spectrum due to non- 
linear evolution as a function of redshift. The non- linear 
Rees-Sciama signal has opposite sign to the ISW on small 
scales where structures are growing more rapidly than in 
linear theory 

For the purpose of estimating the size of effects close 
to the recombination surface I define the lensing poten- 
tial using a visibility source weighting. For the late-time 
contribution this gives results very close to those from us- 
ing a single lens-plane approximation since the thickness 
of last-scattering is very small compared to its distance 
away from us. Using the smooth visibility weighting al- 
lows the importance of near and through-recombination 
contributions to be assessed, though the result is no 
longer directly related to what is observed and the non- 
lcnsing effects discussed in the previous section may also 
become comparable. 

Numerical results are shown in Fig. ^, demonstrating 
that even on large scales it is important to include both 
early and late ISW effects to avoid ler biases, as well 
as large-scale modes and a small correction from the 
Doppler signal. The importance of Rees-Sciama is sig- 
nificantly smaller than claimed in Refs. |$2| [ioj and 
is always a small correction that is safe to ignore in the 
region where the signal is well understood 6 . 



IV. POLARIZATION CORRELATION AND 
BISPECTRUM 



+ 



Ek , lensing 

E(a 




FIG. 5: Correlation of the large-scale .E-polarization with 
the Weyl lensing convergence (which determines the shear), 
the Ricci focussing source —Co, and the sum which deter- 
mines the total magnification correlation. The large-scale E 
polarization is dominated by reionization scattering (optical 
depth r = 0.09 here), and the correlation with the curvature 
perturbation at recombination is small. The shaded band 
shows an indicator of ±1<7 sensitivity to systematic errors 
relative to the cosmic variance limit from lensing, given by 
±y/Cf E C^ + (Cf K ) 2 /l, which is much smaller that the ac- 
tual sensitivity with Planck due to significant large-scale E 
polarization and lensing reconstruction noise. 



bispectrum |T^, g3| that for brevity I shall not reproduce 
here. Anisotropic redshifting does not affect polarization 
components of the bispectra, giving 



b ijk 



c x iT 
< 1 



(s kT cf/ T + 



(19) 



The extension to bispectra involving CMB polarization 
modes is straightforward, giving a Ricci focussing term 
of exactly the same form as Eq. (pil): 



(1*2*3 



h 2dlnx* 



h 



1 3 



(18) 



where Xi , Xj , Xk are either T or E 



There are also 

(larger) standard terms of the Weyl lensing polarization 



The disagreement arises even at the level of the linear power 
spectrum, where Ref. pOj underestimate the lensing correlation 
with linear late-ISW power at / > 100 hy a factor of O(20). As 
this paper was being finalized Ref. ]4l| appeared with results 
that agree with those here. 



The very large-scale polarization comes mostly from 
reionization, and is correlated with the lensing poten- 
tial at up to 30% which gives a Weyl lensing bispec- 
trum contribution that is non-negligible at l\ < 40 [ [l9| . 
There is also a Ricci focussing bispectrum from corre- 
lation of the curvature perturbations at recombination: 
see Fig. ^[ Using polarization the lensing potential can 
be reconstructed with much higher signal to noise than 
with just temperature data, giving constraints on the 
cross-correlation with the temperature and polarization 
that are nearly signal-variance limited and a detection 
at about 8a. Neglecting Ricci focussing to Z max = 2000 
could then bias the total lensing bispectrum amplitude 
by - 0.5a. 

Large-scale polarization is especially useful around 
l\ ~ 60 where the /nl signal in the temperature goes to 
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zero, and for Planck both contributions are small com- 
pared to errors there. It is harder to be sure of the bias 
on /nl from higher l\ without having a full calculation of 
recombination level effects, since a more substantial por- 
tion of the /nl signal comes from li 3> 100, especially 
in the noise-free case. However if we assume that Ricci 
focussing and lensing are only important at l\ < 300, for 
noise free T and E polarization data the Ricci focussing 
bispcctrum gives a bias of /nl ~ 1, corresponding to a 
0.5cr bias on /nl m this case (oy NL ~ 2). See Table. | 
for a summary. The lensing and Ricci focussing signal 
between 100 < l\ < 300 only changes the /nl bias by 
~ 0.2, but dynamical and recombination surface effects 
from smaller scales could conceivably contribute more. 



V. CONCLUSIONS 

Even in the absence of primordial non-Gaussianity, the 
presence of perturbations on a wide range of scales in- 
evitably leads to a non-zero bispectrum in the CMB. By 
directly relating physical scales to observed angles, the 
squeezed limit can be robustly estimated using combi- 
nations of linear-theory results, giving results that agree 
with Refs. Q || in the appropriate limit. The signal is 
dominated by lensing though other terms can give small 
biases to lensing results at the 0.2cr level if ignored, up to 
0.5a with low-noise polarization data. I argued that the 
full-sky approximations presented here should be reliable 
compared to errors for modulation scales l\ < 100. This 
is sufficient for parameter estimation from the lensing 
signal, and local primordial local non-Gaussianity esti- 
mators can be filtered to include only the signal at low 
h 60 with little loss of information for Planck. 

The largest contribution to the bispectrum is from 
lcnsing-tcmpcraturc correlations, which include impor- 
tant contributions from early as well as late ISW. On 
small scales the ISW signal becomes small, and a fully 
self-consistent calculation of all relevant recombination 
effects would be required to calculate the bispectrum. 
It remains possible that at l\ 3> 100 recombination- 
level correlations (e.g. between Doppler temperature 
anisotropics and lenses just in front of last scattering) 
generate detectable net squeezed non-Gaussianities, and 
this should be the subject of future numerical work. The 
same applies for other shapes of non-Gaussianity; only 
squeezed shapes have significant signal to noise involving 
a super-horizon mode, so unfortunately calculating con- 
tributions to other shapes requires detailed study of the 
sub-horizon regime. 

Corrections to the lensing trispectrum (the lensing po- 
tential power spectrum) should be small compared to cos- 
mic variance since the late-time lensing signal is much 
larger. However if the reconstruction is for B-mode lens- 
ing cleaning with high signal to noise, so that the 
actual realization of the reconstructed lensing potential 
matters, small errors due to Ricci focussing and dynam- 
ical anisotropics at recombination may require more de- 



tailed study. 
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Appendix A: Linear temperature anisotropies from 
anisotropic redshifting 

In the delta-function visibility approximation the 
cqual-tcmpcrature recombination surface is observed by 
us to have anisotropics due to differing amounts of red- 
shift in different directions. The result is very well known 
(for a more detailed review see e.g [^3) ) , but I recap here 
for completeness and rephrase in a somewhat unfamil- 
iar way to demonstrate the close analogy with the ray 
bundle focussing result. 

In the conformal Newtonian Gauge, integrating the 
photon geodesic equation gives the energy today E A in 
terms at the energy at a (conformal) time n 



a A E A = a{rf)E{rf) 



1 + *(?/) - * y 



(Al) 

where energies are those observed by observers with no 
(Newtonian gauge) peculiar velocity and the ISW inte- 
gral is along the line of sight. The potentials at -q and A 
appear local, but arise from integrating total derivative 
terms along the photon path. This can equivalently be 
written 



a A E A = a{rj)E{rf) 



1 - $(?/) + 



>IA 



dr]dx(^ + $) 



(A2) 

where x is the comoving radial distance and now the 
integral term involves spatial gradients (and hence is not 
zero in matter domination). For a constant temperature 
surface the time of emission is rj = ?/„ + 5r), with Si] = 
—Spy/ p 7 = A 7 /(4W), so the scale factor at the perturbed 
emission time is a(rj) = a*(l + A 7 /4). Accounting for 
Doppler shifts due to peculiar velocities and converting 
to the fractional temperature perturbation then gives the 
temperature anisotropy 

A f riA 
AT(n) = — y - - $ + n • (v A - v) + / drj n ■ V(* + $) 



C 7 + n • (v A - v) + / dr/n • V(* + 
Jn 



$) (A3) 



where I defined a A (l — $^) = 1 and C 7 is the gauge- 
invariant 3-curvature on the constant-temperature re- 
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combination hypcrsurfacc. On large scales there is a 
cancellation between £ 7 which is positive in an overden- 
sity, describing getting relatively hotter due to the Ricci 
defocussing leaving an overdensity (relative space con- 
traction, so gets a relative blue shift), and the effect of 
climbing out the potential well: a line-of-sight integral 
over the line-of-sight gradients in the potentials. The 
net effect in matter domination is the Sachs- Wolfe result 
$/3. 



The form of this equation is very similar to the re- 
sult for the beam area, Eq. (|J), involving the curvature 
perturbation at the source, and an integral of potential 
gradients (the radial analogue of the transverse gradients 
that give rise to convergence and shear). That they both 
depend on £ 7 should be unsurprising since photon energy 
and beam size are both affected by the expansion of space 
(C 7 = 5N = 8 In a). 
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